Finding out the different types of land use layers in the landcover column
unique(hawai_land$landcover) #Informs what types of land use covers exist in the dataset
#Creating a new column in the data set and categorizing the following land use covers in order to streamline organization
landmap <- hawai_land %>%
mutate(
landmapnew = case_when(
landcover == "Cropland and Pasture" ~ "Farmland",
landcover == "Commerical and Services" ~ "Urban",
landcover == "Residential" ~ "Human Dwellings",
landcover == "Evergreen Forest Land" ~ "Forest",
landcover == "Other Urban or Built-up Land" ~ "Urban",
landcover == "Mixed Rangeland" ~ "Rural",
landcover == "Industrial" ~ "Urban",
landcover == "Streams and Canals" ~ "Watershed",
landcover == "Orchards, Groves, Vineyards, Nurseries and Ornamental Horticultural Areas" ~ "Rural",
landcover == "Shrub and Brush Rangeland" ~ "Rangeland",
landcover == "Forested Wetland" ~ "Forest",
landcover == "Reservoirs" ~ "Watershed",
landcover == "Nonforested Wetland" ~ "Forest",
landcover == "Bare Exposed Rock" ~ "Rock",
landcover == "Sandy Areas Other than Beaches" ~ "Sand",
landcover == "Transportation, Communications and Utilities" ~ "Urban",
landcover == "Herbaceous Rangeland" ~ "Rangeland",
landcover == "Beaches" ~ "Beaches",
landcover == "Other Agricultural Land" ~ "Farmland",
landcover == "Lakes" ~ "Watershed",
landcover == "Strip Mines, Quarries, and Gravel Pits" ~ "Mining",
landcover == "Mixed Barren Land" ~ "Barren Land",
landcover == "Bays and Estuaries" ~ "Watershed",
landcover == "Mixed Urban or Built-up Land" ~ "Urban",
landcover == "Transitional Areas" ~ "Unknown",
landcover == "Industrial and Commercial Complexes" ~ "Urban",
landcover == "Confined Feeding Operations" ~ "Urban",
landcover == "0" ~ "NA"
)
)
tmap_mode("view")
Creating the map which contains all the new categories for the land use layers in Hawaii
tm_shape(landmap) +
tm_polygons("landmapnew", title = "Land Cover Types for All Main Hawaiian Islands", style = "fixed") +
tm_borders() +
tm_layout(legend.title.size = 1,
legend.text.size = 0.6,
legend.position = c("left","bottom"),
legend.bg.color = "white",
legend.bg.alpha = 1) +
tm_basemap("Esri.NatGeoWorldMap")